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CHAPTER  I .  INTRODUCTION 


1 .1  Background 

The  resistive  boundary  condition  was  devised  In  order  to 
simplify  the  solution  of  boundary  value  problems  Involving  thin  lossy 
dielectrics,  and  Harrington  and  Mautz  (1975)  have  shown  that  It 
provides  a  good  model  for  such  materials.  Practical  examples  of 
"resistive  materlils"  range  from  man-made  lossy  composite  materials 
used,  for  example,  for  the  fin  of  a  missile  to  reduce  Its  radar 
cross  section,  to  snowflakes  and  smoke  particles.  A  precise  definition 
of  the  resistive  boundary  condition  will  be  given  later  In  this 
chapter. 

A  number  of  dielectric  structures  have  been  simulated  by  this 
boundary  condition.  Infinite  structures  such  as  the  half  plane,  strip 
(Senior,  1979a, b)  and  wedge  (LaHale,  1981)  have  been  treated 
analytically,  using  transform  techniques.  Closed  shells  (Harrington 
and  Mautz,  1975;  Senior  and  Well,  1977)  were  treated  by  numerical 
methods.  There  have  been  no  previous  numerical  analyses  of  finite, 
open  resistive  structures  because  of  the  Inherent  numerical  difficulty 
associated  with  such  structures.  This  difficulty  Is  common  to  any 
finite,  open,  thin  structure,  resistive  or  perfectly  conducting 
alike.  It  arises  because  only  the  highly  singular  E-fleld  Integral 
equation  Is  applicable  and  the  situation  has  been  aggravated  by  the 
common  belief  that  only  E-fleld  formulations  derivable  from  the 
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Franz  representation  are  valid  (Stratton,  1941;  Bouwkamp,  1947).  The 
various  representations  and  their  numerical  Implications  will  be 
discussed  In  Chapter  II.  Because  the  source  of  numerical  difficulty 
In  the  treatment  of  plates  Is  the  same  for  all  boundary  conditions. 

It  Is  helpful  to  review  the  techniques  used  for  treating  perfectly 
conducting  plates. 

Mlttra  et  al .  (1973)  and  Rahraat-Saraml  and  Mlttra  (1974)  suggested 
a  way  of  partially  Integrating  an  E- field  equation,  thus  reducing  the 
order  of  Its  singularity.  Their  method  Is  a  two-dimensional  analog 
of  Hallen's  Integral  equation  for  wire  antennas.  Unlike  the  one¬ 
dimensional  case,  the  general  homogeneous  solution  of  the  differential 
operator  Is  not  known.  Mlttra  et  al .  proceeded  by  assuming  that  the 
unknown  general  solution  Is  expandable  In  a  Fourier  series  and 
determined  the  coefficients  of  the  expansion  by  fmposing  an  edge 
condition  along  the  rim  of  the  plate.  Their  method  works  for  perfectly 
conducting  plates,  but  does  not  extend  to  resistive  ones. 

Recently,  Glisson  and  Wilton  (1980)  came  up  with  a  different 
method.  They  used,  but  without  justification,  an  E-fleld  Integral 
equation  based  on  the  Stratton-Chu  representation.  As  will  be  seen 
In  Chapter  III,  their  numerical  scheme  guarantees  both  continuity 
of  the  vector  potential  and  satisfaction  of  the  edge  condition.  Their  . 
method  Is  a  breakthrough  from  a  numerical  viewpoint,  but  requires  a 
sounder  theoretical  basis. 

The  Bablnet  complement  of  a  perfectly  conducting  plate  is  an 
aperture  In  a  metallic  screen.  Because  of  its  great  Importance,  the 
problem  of  diffraction  by  an  aperture  In  a  plane  screen  has  been 
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Investlgated  by  many  authors.  Most  treatments  Involve  some  kind  of 
an  approximation,  either  physical  optics  (cf.  Bouwkamp,  1953),  the 
geometrical  theory  of  diffraction  (Keller,  1962)  or  low  frequency 
expansions.  The  last  method  Is  an  example  of  Rayleigh  scattering 
(Lord  Rayleigh,  1897).  A  systematic  procedure  for  treating  diffraction 
by  small  obstacles  was  developed  by  Stevenson  (1953).  In  the  many 
papers  that  followed,  a  variety  of  low  frequency  diffraction  problems 
were  solved.  Van  Bladel  (1977)  provides  a  comprehensive  list  of  the 
work  published  through  1975. 

Integral  equations  for  the  fields  of  small  apertures  were 
developed  by  Bethe  (1944).  Van  Bladel  (1970)  and  his  colleagues 
(De  Neulenaere  and  Van  Bladel,  1977;  De  Smedt  and  Van  Bladel,  1980) 
attempted  to  solve  them  numerically,  but  with  only  limited  success. 
Butler  (1974)  Integrated  those  Integral  equations  In  a  manner 
analogous  to  that  used  by  Mlttra,  and  numerical  results  have  been 
reported  by  Umashankar  and  Butler  (1974).  A  different  approach  based 
on  a  finite  element  formulation  was  used  by  Okon  and  Harrington  (1980). 

1 .2  Outline  of  the  Work 

The  task  Is  to  determine  the  fields  scattered  by  a  resistive 
plate  S,  lying  In  the  x-y  plane  (Fig.  1.1),  when  illuminated  by  a 
plane  wave: 

jlnc  .  *  e1kk-r  .  jlitt  .  Yietkk-r  (,.i) 


l 


-5- 


The  time  convention  e 


-lut 


Is  assumed  throughout. 


The  unit 


vector 


*  *  4  4  A  4  <1  A  4 

k  ■  -(x  cos  4  cos  4  +  y  cos  t  sin  $  +  z  sin  <fr  )  (1.1) 

Is  In  the  direction  of  propagation  of  the  Incident  wave,  a,  b  and  k 
form  a  right-handed  orthonorraal  set.  The  resistivity  R  of  the  plate 
Is  In  general  complex  and  a  function  of  position.  Its  value  can  be 
allowed  to  go  to  zero  (the  perfectly  conducting  limit). 

Th*  problem  Is  solved  numerically  by  applying  the  moment 
met'cJ  to  a  suitable  Integral  equation  for  the  induced  currents. 
Because  of  the  finite  capacity  of  any  computer,  the  largest  plate 
that  can  be  treated  at  present  is  about  a  square  wavelength  in  area. 
Although  this  Is  not  large  enough  to  permit  direct  comparison  of  the 
results  with  those  obtained  by  asymptotic  methods  such  as  the 
geometrical  theory  of  diffraction,  plates  of  this  size  are  of 
practical  Interest  and  display  similar  characteristics  to  larger 
plates. 

Special  attention  is  paid  to  the  low  frequency  limit,  where 
semi-analytic  results  can  be  derived  from  the  Rayleigh  series.  The 
first  two  terras  of  the  expansion  are  determined,  and  efficient  methods 
are  developed  for  the  solution  of  the  associated  integral  equations. 

The  resistive  boundary  condition  and  its  relation  to  other 
types  of  boundary  conditions  are  discussed  in  the  remaining  section 
of  this  chapter.  The  first  half  of  Chapter  II  Is  devoted  to  a  dis¬ 
cussion  of  the  representation  of  an  electromagnetic  field.  A  proof 
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of  the  adequacy  of  the  Stratton-Chu  formulation  for  a  plate  is 
provided  there.  The  method  of  moments  and  some  of  Its  ramifications 
are  then  presented,  and  these  two  topics  lay  the  foundation  for 
the  numerical  solution  of  the  problem  In  Chapter  III.  Several 
Integral  equations  are  derived  and  evaluated.  One  is  selected  for 
numerical  treatment  on  the  basis  that  it  provides  the  fastest 
convergence.  Results  are  compared  with  experimental  data  and  their 
Implications  are  discussed. 

Chapters  IV  and  V  deal  with  the  low  frequency  limit.  The 
analytical  derivations  are  carried  out  In  Chapter  IV  and  the 
numerical  solution  In  Chapter  V.  Comparison  with  other  methods  and 
some  experimental  and  theoretical  bounds  for  the  dipole  moments  are 
also  covered  In  the  latter  chapter. 

Concluding  remarks  and  suggestions  for  the  extension  of  this 
work  are  reserved  for  the  last  chapter. 

1 .3  The  Resistive  Boundary  Condition 

A  resistive  sheet  is  characterized  by  three  properties.  It  is 
infinitesimally  thin,  carries  only  electric  currents  and  these  are 
proportional  to  the  tangential  components  of  the  electric  field. 

Mathematical  expressions  for  these  conditions  are: 

0*  =  nj|*  =  0  (1.3) 

and 

n  x  E  *  Rn.3  *  Rn„(n^H)|*  .  (1.4) 


The  proportionality  factor  R  Is  related  to  the  material 


properties  by 


R  -  11m  ^  »  11m  fleJ-JT)F  0-5) 

a  -*■  »  er  •  ** 

h  *►  o  h  -+•  o 

where  h  Is  the  thickness  and  a  the  conductivity  of  the  material. 

In  practice,  this  boundary  condition  Is  applied  to  nonmagnetic 
(y  3  uQ),  highly  conductive  materials  whose  thickness  h  Is  very  small 
compared  with  the  wavelength  x  and  the  penetration  depth  <5.  The 
material  Is  thus  partially  transparent.  The  value  R  *  0  corresponds 
to  perfect  conductivity. 

Equation  (1.5)  does  not  determine  the  frequency  dependence  of 
R.  Empirically,  the  resistivities  of  resistive  sheets  remain 
constant  over  a  broad  range  of  frequencies,  and  In  fact  are  measured 
at  dc.  Thus  It  can  be  assumed  that  the  resistivity  R  is  Independent 
of  k,  at  low  frequencies  at  least. 

The  resistive  boundary  condition  is  related  to  two  other  types 
of  boundary  conditions:  the  conductive  and  the  impedance  boundary 
conditions. 

The  conductive  boundary  condition  is  the  dual  of  the 
resistive  one.  A  conductive  sheet  supports  only  magnetic  currents 
(J  *  n,H|_  *  0),  which  are  related  to  the  tangential  components  of 
the  magnetic  field  by 
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fLH  -  R*n*5*  ■  -R*n.(n„E)|* 

and 

R*  \  • 
n  -*•  o 

This  Is  a  fictitious  kind  of  material.  The  only  practical  importance 
of  the  boundary  condition  Is  that  combined  with  a  resistive  sheet 
It  produces  an  Impedance  sheet. 

An  impedance  sheet  Is  an  opaque  (l.e.,  h  »  &),  thin  (h  <<  a) 
layer  supporting  both  electric  and  magnetic  currents.  These  currents 
are  related  to  each  other,  and  on  each  face  of  the  sheet 

n.Cn.E)  *  -nZn„H  .  (1.6) 

For  a  planar  sheet  this  relation  reads 

Ex  *  ^ZHy  *  EJ  “  ±nZHx  ‘ 

The  upper  and  lower  signs  refer  to  the  top  and  bottom  faces 
respectively. 

By  writing  the  resistive  and  conductive  boundary  conditions 
for  planar  sheets  in  the  forms 

E*  +  £  -  -2R(H*  -  H~)  ,  E*  +  E~  =  2R(H*  -  H‘) 

c  -  ■  -m  *  h;>  •  £y  -  e;  •  h*  <hx  ♦  «;> 
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and  assuming  coexistence  of  magnetic  and  electric  currents,  Senior 
(1978)  obtained 


These  relations  turn  Into  the  impedance  boundary  condition  If 
nZ  =  2R  =  1/2R*. 

Thus  the  electric  and  magnetic  currents  Induced  on  an  Impedance 
sheet  are  the  same  as  those  on  a  resistive  and  a  conductive  sheet 
with  R  *  (l/2)nZ  and  R*  *  l/(2nZ),  respectively. 

Having  solved  for  the  currents  on  a  resistive  plate,  the 
magnetic  currents  on  a  corresponding  conductive  plate  can  be 
obtained  by  duality.  The  Impedance  plate  Is  then  the  superposition 
of  both  types  of  current  sheets. 
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CHAPTER  II.  THEORETICAL  BACKGROUND 


2.1  Representations  of  the  Electromagnetic  Field 

Integral  representations  are  powerful  tools  in  the  theory  of 
electromagnetism.  They  provide  a  means  for  the  derivation  of 
integral  equations  and  of  asymptotic  low-  and  high-frequency 
approximations. 

Several  representations  are  available.  Although  they  will 
be  shown  to  be  equivalent,  each  possesses  certain  properties  that 
make  it  advantageous  for  some  applications  and  not  for  others. 

A  couple  of  cases  will  demonstrate  this  point. 

The  total  field  is  the  sum  of  the  scattered  and  incident 
fields.  Since  the  latter  is  given,  only  the  scattered  field  will 
be  considered. 

The  most  widely  used  representation  is  due  to  Stratton  and 
Chu  (Stratton,  1941;  p.  464).  It  expresses  the  scattered  field 
in  terms  of  the  surface  values  of  both  the  electric  and  magnetic 
fields  on  the  surface  S  of  a  scattering  body: 


is  =  ikZ  f  n'.H  G  dS ’  -  v  f  n'-E  G  dS ’  +  v.  f  n'.E  G  dS '  (2.1) 

Js  Js  ^S 


and 


H*  =  -ikY  j  n'.E  G  dS '  -  vf  n'-H  G  dS '  +  7.  f  n'.H  G  dS '  .  (2.2) 

Js 


Identifying  n  x  H,  n  x  E,  n*E  and  n*H  with  J,  J*,  p  and  p* 
respectively,  the  various  Integrals  In  (2.1)  and  (2.2)  are  recognized 
as  vector  and  scalar  potentials.  Hence  the  fields  can  be  written  as: 


Es  *  1 kZA  -  7$  +  v^A* 


(2.3) 


H  *  -IkYA*  -  V4>*  +  7. A 


(2.4) 


These  two  formulations  are  entirely  equivalent. 

The  charges  and  currents  are  related  and  the  relations  are 
given  by  the  continuity  equations: 


ikY  £■ 


(2.5) 


-ikZ 

y 


(2.6) 


These  equations  supplement  (2.1)  and  (2.2),  or  (2.3)  and  (2.4)  to 
form  the  complete  representation. 

Another  formulation  attributed  to  Franz  (Jones,  1964)  requires 
knowledge  of  the  induced  currents  alone: 


Es  =  (k2  +  V7-)  j  n'  H  G  dS'  +  7  f  n' J  G  dS ’  ,  (2.7) 


P-  (k2  +  77-)  f  n’  E  G  dS'  +  7  f  n\ 


H  G  dS*  .  (2.8 
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The  terras  of  this  representation  are  electric  and  magnetic 
Hertz  potentials  (linearly  related  to  the  vector  potentials).  An 
unconventional  derivation  of  the  representation  which  avoids  the 
use  of  Green's  theorem  Is  presented  In  Appendix  A. 

It  Is  well  known  that  the  fields  represented  by  the  Franz 
formulation  are  Maxwellian.  Stratton-Chu  fields,  however,  satisfy 
Maxwell's  equations  if  and  only  if  they  are  equivalent  to  those  of 


Franz. 


By  applying  Stokes'  theorem  and  Maxwell's  equation 


V^H  =  -ikYE  to  /  V'„[HG]  -n  dS '  it  is  found  that  (ibid) 


i'.H  G  dS 


'  =  -ikYv  J  n ' -i 


G  dS'  +  d>  7'G  H-dc 

Jc 


(2.9) 


An  analogous  result  holds  for  the  dual  case  (E H,  H  -E,  Y  -«■  Z). 
It  is  thus  evident  that  the  two  representations  are  equivalent  when 
the  line  integrals  vanish.  The  task  is  therefore  to  show  that  the 
line  integrals  a* e  indeed  zero. 

The  case  of  a  smooth  closed  surface  S  is  trivial.  The  closed 
path  C  degenerates  into  an  arbitrary  single  point  of  S  and  the 
line  integrals  disappear. 

Investigation  of  the  behavior  of  the  fields  in  the  vicinity 
of  an  edge  establishes  the  same  result  for  surfaces  with  edges. 
Attaching  a  local  cylindrical  coordinate  system  p,$,c  (Fig.  2.1)  and 
imposing  the  condition  that  the  energy  stored  in  the  neighborhood 
of  the  edge  is  finite,  the  fields  are  seen  to  be  0(pq_1)  with  q  >  0. 


I 


The  p  and  41  components  of  Maxwell's  equation  7  x  H  ■  -IkYE 


3H  3H 


p  3$  3C 


(2.10) 


3H  3H 


at  ap 


-ikYE, 


(2.11) 


If  all  components  of  E  and  H  are  0(pq"  ),  then  the  first  term  on  the 
left  hand  side  of  (2.10)  and  the  second  term  on  the  left  hand  side  of 
(2.11)  are  0(pq~2).  Since  the  other  terms  are  0(pq_1)  It  follows  that 


1 

11m  -  77^ 


p  -*■  0 


p  34 


(2.12) 


=  0  . 


(2.13) 


p  0 


The  other  possibility  Is  that  H_  Is  O(p^)  In  which  case  11m  Hz  «  0. 

p  -►  0 

Equations  (2.12)  and  (2.13)  have  two  distinct  solutions:  one 
has  q  t  1 ,  implying  lim  Hz  =  0  (the  other  components  can  be  singular), 

p  -*■  0 

and  the  other  has  Hz  independent  of  both  p  and  4 ,  implying  q  3  1. 

In  both  cases,  Hz  is  continuous  at  the  edge,  Independently  of  the 

boundary  condition  that  holds  on  S.  In  particular,  if  S  3  S  U  S 

1  2 

where  S  and  S  are  (open)  smooth  surfaces  that  meet  at  an  angle 
1  2 

along  C 
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J  7'.(HG)  •  n  dS'  «  (£(ff  -  ff  )G  •  dc  (2.14) 

s  u  s  c  1  2 

1  2  U 

holds.  IT  and  W  are  the  values  of  H  on  S  and  S  respectively. 
12  12 

Since  H  Is  Independent  of  4  the  Integral  on  the  right-hand  side 

M 

of  (1.14)  vanishes.  This  result  can  be  generalized  when  S  *  U  S 

m*l  m 

and  S  has  several  edges. 

The  continuity  of  H  also  guarantees  the  disappearance  of  the 

line  Integrals  In  the  case  of  an  aperture  In  a  perfectly  conducting 

screen.  The  Integral  /v'GH  *di;  vanishes  on  the  screen,  hence  it 

C 

vanishes  on  the  aperture  side  of  the  edge  also. 

In  the  next  section  It  will  be  shown  that  the  choice  of 
representation  has  a  profound  effect  on  the  convergence  properties 
of  a  numerical  solution. 

2.2  Numerical  Solution  of  Integral  and  Differential  Equations  by 
the  Moment  Method 

2.2.1  General  Procedure.  The  method  of  moments  (cf.  Harring¬ 
ton,  1968)  was  devised  to  solve  numerically  equations  of  the  form: 

Lf  *  g  ;  in  n  C  RP  ,  (2.15) 

where  L  Is  an  operator  mapping  the  unknown  function  f  €  X  (over  R^) 
Into  the  known  function  g  €Y  (over  R^). 

An  Inner  product  Is  defined  In  Y  as 

<u,v>  *  J  uv  dx^  ;  )/u,v  f  Y  (2.16) 

n 


H 

A  set  {w  }  C  Y  of  weighting  (or  sampling)  functions  Is  selected. 
Taking  the  Inner  product  of  the  given  equation  with  wm  we  get 

<w_,Lf>  ■  <w  .g>  ;  1  <  m  <  M 

m  m  —  — 

The  function  f  is  then  expanded  In  terras  of  a  basis  {f  }“  of  X 

n  1 

oo 

f  ■  I  Vn  <2-17> 

n*i 

and  the  series  is  truncated  at  some  finite  N.  The  result  Is  a  set 
of  M  algebraic  equations  for  the  N  unknown  coefficients: 

N 

£  <wm’LV  an 

n=i 

or,  in  matrix  notation, 

!  a  .  (2.i9) 

where  Wmn  =  <wm,Lfn>  and  ^  =  <wm  ,g>.  The  system  of  equations 
(with  M  =  N)  is  solved  by  Gaussian  elimination  and  backwards 
substitution. 

2.2.2  Subsectional  Bases  and  C°[n]  Continuity.  The  choice 
of  the  basis  {fnJ  in  which  f  is  expanded  is  a  crucial  step  in  the 
solution.  Obviously  a  set  of  basis  functions  which  resemble  the 
required  solution,  that  is,  satisfy  the  same  boundary  conditions  and 
have  the  same  continuity  properties,  will  lead  to  rapid  convergence. 
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and  only  a  few  terras  of  the  expansion  are  then  needed  to 
describe  f. 

Generally  It  Is  difficult  to  "guess"  a  good  basis,  and  It  Is  a 
almost  Impossible  task  for  a  domain  n  of  complex  geometry  with 
complex  boundary  conditions.  It  Is  much  easier  to  divide  G  Into 
subdomains  such  that 

N 

fl  ■  U  nN  (2.20) 

n=i 

and  use  basis  functions  which  have  support  on  nn: 


Such  a  basis  is  called  a  subsectional  basis. 

Quite  often  the  nature  of  the  operator  L  dictates  that  f 
possess  certain  continuity  properties.  A  subsectional  basis  may  lead 
to  a  representation  of  f  which  does  not  possess  the  required 
continuity,  and  thus  results  in  poor  convergence. 

For  a  planar  domain  n  C  R2  there  is  more  than  one  choice  of 
a  subsectional  basis  that  guarantees  C°[n]  continuity.  One 
possibility,  which  will  be  extensively  used  in  this  work,  uses 
rectangular  cells: 
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*t  •  /„  -  f  i  y  i  *„  ♦  f] 


(2.22) 


and  basis  functions  that  are  the  tensor  product  of  two  quadratic 
polynomials  In  s  and  n  respectively,  c  =  x  -  xm,  n  =  y  -  yfl: 


a  +oC+an+a5n+aC2 
1  2  3  4  5 

+  a  n2  +  a  £zn  +  a  £n2  +  a  5zn 
6  7  8  9 


;2n2  ;  (x.y)  €*  nB 

;  (x.y)  f  nn 


(2.23) 


The  term  In  £2n2  may  eventually  be  dropped,  but  it  is  more 

convenient  to  retain  it  for  the  time  being. 

The  coefficients  (which  depend  on  m,n)  are  selected  so  that 

f  attains  its  exact  value  at  nine  predetermined  points  of  Qmn-  It 

is  shown  in  Appendix  B  that  If  these  points  are  taken  to  be  the 

four  corners  of  £J  the  four  midpoints  of  its  sides  and  Its 
mn 

center,  then 


f  *  I 


With  this  choice  of  nodes  In  flmn  (Fig.  2.2a) 


it 

3X 


a 

2 


flxm*  r 


,yn*  "  f^xm  '  2*  ^ 
Ax 


ill 

ay  I 


f(Vyn  +  Z~  }  -  f(Vyn  '  2  -  ) 


(Vyn} 


Ly 


f(*m  +  r»  yn}  *  f(xm  ~  ?'  yn}  ~  2f(V 

AX2 


AX 


a2f 

ax2 


^xm’yn^ 


2a  *  4 

5 


a2f 

ay2 


2a 


f{Vyn  +  Z~  ]  +  f(Vyr 


-  ^  ) 
2  ; 


2«V 


Ay2 


(2 

and  the  right-hand  sides  are  finite  difference  representations  of 

these  derivatives,  numerically  equivalent  to  (2.17)  with  (2.23). 

The  expressions  given  by  (2.24)  are  particularly  useful  when  f  is 

a  weighted  Integral  of  another  unknown  function  to  be  determined. 

The  choice  of  sampling  points  is  crucial.  The  dramatic 

effect  of  loss  of  continuity  as  the  points  1,2,3,  and  4  are  moved 

away  from  the  boundary  of  n  will  be  demonstrated  when  numerical 

mn 


results  are  discussed. 


Fig.  2.2  (a)  The  Nodes  of  a  Rectangul 

General  Quadrilateral  Cell. 
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Before  concluding  this  section.  It  Is  appropriate  to  remark 
that  the  rectangular  cell  can  be  mapped  Into  a  general 
quadrilateral  cell.  A  quadratic  transformation,  based  on  the  same 
nine  points.  Is  a  very  convenient  means  to  effect  the  mapping  at 
least  approximately  (Fig.  2.2b).  The  applicability  of  the  method 
Is  thus  extended  to  more  general  geometries. 

2.2.3  Subsectional  Bases  and  C*[n]  Continuity.  Some 
operators,  such  as  32/3x3y,  Impose  higher  continuity  constraints 
on  f,  and  this  particular  operator  requires  C‘[fl]  continuity.  It 
is  much  more  difficult  to  achieve  this  kind  of  continuity  with 
subsectional  bases.  Both  f  and  Its  partial  derivatives  must  be 
specified  at  the  nodes  of  n  .  Even  then,  no  simple  polynomial 
expression  which  ensures  continuity  of  f  and  its  first  derivatives 
Is  possible  (Zlenklewicz,  1977,  Chapter  10). 

Many  attempts  have  been  made  to  overcome  this  difficulty, 
with  different  degrees  of  success  [ibid.  Chapters  10-12].  For 
mechanical  problems  it  has  been  recognized  that  the  best  approach 
is  to  reformulate  the  problem  so  that  only  C°[n]  continuity  Is 
required.  The  original  equation  (2.15)  is  replaced  by  a  set  of 
equivalent  equations  involving  operators  of  lower  order.  It 
turns  out  that  It  Is  more  economical  to  Increase  the  number  of 
equations  than  to  use  more  parameters  per  cell. 


i 
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CHAPTER  III.  SCATTERING  BY  RESISTIVE  PLATE-DYNAMIC  FORMULATION 


3.1  Integral  Equations 

Integral  equations  for  the  induced  currents  are  derived  from 
the  field  representations  by  taking  the  observation  point  to  He  on 
the  surface  S  of  the  scatterer  and  applying  the  boundary  conditions 
there. 

Several  E-field  formulations  are  possible  but  the  H-field 
equations  becomes  an  identity  for  an  infinitely  thin  open  body. 
Selection  of  the  preferred  equation  for  numerical  treatment  is  based 
on  the  conclusions  of  the  preceding  chapter.  From  a  numerical 
standpoint,  the  equation  that  imposes  the  lowest  order  of  continuity 
requirements  is  the  best. 

One  equation  is  deduced  from  the  Franz  representation  using 
the  boundary  conditions  (1.3)  and  (1.4): 


Rz.j  =  z„E1nc  +  ^  zjk2  +  vv*) 


1 


J  G  dS 1 


(3.1) 


The  dyadic  operator  vv  contains  second  order  mixed  derivatives  such 
as  32/3x3y  and  requires  C 1 [S J  continuity  of  the  vector  potential. 

This  is  hardly  encouraging  and  the  need  for  other  formulations  is 
evident. 

Operating  on  one  Cartesian  component  of  (3.1)  with  32/3x2  +  k2 


and  with  32/3y2  +  k2  on  the  other,  subtracting  and  simplifying 
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the  incident  field  term  using  all  three  relations 

v  •  Einc  =  0  ,  v  x  E1nc  =  1kZH1nc 

and 

(v2  +  k2)E1nc  =  0 


leads  to  an  equation  for  one  component  of  3.  Repeating  the  whole 
process  with  the  two  differential  operators  Interchanged  results  In 
an  equation  for  the  second  component.  Combined  together  they 
provide  the  Integral  equation: 


ikZ(v2 


/ 


+  k2)  J  G  dS'  =  1kZz„ 


3H 


inc 


3Z 


+  R[k2  +  z.7(z-7j]J  .  (3.2) 


This  is  valid  only  if  R  Is  constant,  but  the  equation  has  even  more 
serious  shortcomings.  C1  continuity  of  3  itself  is  now  required. 

■  ^  n/* 

The  equation  fails  at  grazing  incidence  since  3H  /3z  introduces 
a  factor  sine^  into  the  excitation  term.  The  transition  to  a 
perfectly  conducting  plate  is  also  discontinuous  and  when  R  ■  0 
there  is  no  coupling  between  the  two  components  of  the  currents. 
Satisfaction  of  the  boundary  condition  z*H  »  0  (on  S)  is  no  longer 
guaranteed  and  must  be  imposed  separately. 

The  only  advantage  this  formulation  may  offer  is  the 
possibility  of  partial  integration.  When  R  ■  0  this  equation  was 
integrated  by  Mittra  et  al . ,  (1973)  by  using  Green's  theorem  in  a 
plane  and  the  two-dimensional  "free  space"  Green's  function.  Their 
method  is  capable  of  removing  the  differential  operator  v|  +  k2. 
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However,  it  is  impossible  to  integrate  z„v(z*v„)  at  the  same  time, 
and  thus,  if  R  t  0  we  are  still  left  with  the  requirement  for  C1 
continuity  of  3. 

A  third  set  of  equations  can  be  derived  from  the  Stratton-Chu 
representation.  The  boundary  conditions  (1.3)  and  (1.4),  when 
substituted  into  (2.1),  produce  the  integral  equation: 

RzJ  =  zj1nc  +  ikZz^Jj  G  dS '  -  zAV  J G  dS '  (3.3) 

subject  to  the  constraint 


(3.4) 


The  highest  order  continuity  required  is  C°[S],  and  the  burden  of 
continuity  is  now  shared  between  the  scalar  potential  (in  3.3)  and 
the  current  (in  3.4).  The  price  paid  is  the  increased  number  of 
equations  that  have  to  be  solved  simultaneously.  In  a  later  section 
it  is  shown  how  that  difficulty  is  circumvented. 

Having  a  valid  and  convenient  formulation,  it  is  feasible  now 
to  discuss  such  matters  as  uniqueness  of  the  solution  and  numerical 
procedures . 

3.2  Uniqueness  of  the  Solution 

Jones  (1974)  has  examined  the  question  of  uniqueness  of  the 


solution  in  the  case  of  a  perfectly  conducting  solid  body.  He 
showed  that  E-field  integral  equations  have  more  then  one  solution 
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at  wave  numbers  which  are  eigenvalues  corresponding  to  the  TE  modes 
of  a  cavity.  Identical  In  size  and  shape  to  the  body  and  having 
perfectly  conducting  walls. 

Regarding  the  plate  as  a  limit  of  a  hollow  right  prism  of 
altitude  h,  and  repeating  Jones'  argument,  we  see  that  if  the 
equations  (3.3)  and  (3.4)  have  more  than  one  solution,  the 
differences  J,p  of  these  solutions  satisfy 

RnJ  =  ikZruJj  G  dS*  -  n.v  G  dS  ‘ 
and 

v .  3  .  f  *  . 

These  are  the  equations  satisfied  by  the  eigenmodes  of  the 
cavity.  All  the  eigenvalues  of  the  equations  are  complex,  but 
even  so  the  accuracy  and  stability  of  the  solutions  for  k  real  may 
be  seriously  affected  by  the  presence  of  a  nearby  complex  pole. 
Every  mode  of  these  equations  can  be  put  into  a  one-to-one 
correspondence  with  a  TE  mode  of  a  perfectly  conducting  cavity 
of  the  same  geometry.  As  the  altitude  h  goes  to  zero  all  TE  modes 
disappear.  Hence  the  solution  of  (3.3)  and  (3.4)  is  unique  for  any 
value  of  R,  including  R  =  0. 
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3.3  Numerical  Solution 

Equations  (3.3)  and  (3.4)  are  solved  by  tbe  moment  method. 
Written  in  scalar  form,  these  equations  are: 


RJx  =  eJ"c  +  ikZ  Jjx G  dS'  -  J  6  dS* 
RJy  =  E™c  +  ikZ /  V  dS'  -  J*-  G  dS‘ 
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If  the  finite  difference  expression  (2.22)  is  to  be  applied  to 

(3.4),  it  is  clear  that  the  sampling  points  of  J  and  J  must  be 

x  y 

shifted  by  ax/2  and  Ay/2,  in  the  x  and  y  directions,  respectively, 
relative  to  the  sampling  point  of  p.  Bearing  in  mind  that  it  is 
inadvisable  to  sample  p  on  the  edges  of  the  plate  (where  p  is 
singular  if  R  *  0),  and  that  Jx  and  are  zero  on  edges  parallel 
to  the  y-  and  x-axes,  respectively,  one  arrives  at  Glisson  and 
Wilton's  (1980)  scheme  of  interwoven  meshes  (Fig.  3.1)  for  a  rectangular 
plate. 

A  corresponding  expansion  of  the  unknown  functions  is 


J 

x 


I 


P=i 


N 


I 


rect 

P  +1/2, q 


x 


X42.+1  /-2.»q 

AX 


(3.5) 


i 


1 


l 


p=i  q=i 


£.  = 


e 


M  N 


I  I 


rect 


rect 


(3.7) 


where 


Jx.  =  J  =  J  =  J  =0 
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from  the  edge  condition.  Using  sampling  functions 
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for  (3.3a),  (3.3b)  and  (3.4),  respectively,  and  using  (2.22), 
equations  (3.3)  and  (3.4)  are  converted  to 
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An  expansion  of  the  potential  corresponding  to  (3.5)  through 
(3.7)  is 


M-i  N 


ra+1 /2,n 


J.  G, 


x“m+l/2,n  ^  X-  2L  ’x. 


^  P+l/2,q 


bm+l/2,n,p+1/2,q  etc. 


where 


y”y 


and 


1  »  V  9  K  |  X  / 


G  dS  * 

y  ,v 


K  »X 


Using  the  translational  symmetry  of  G: 


m+a,n+8,p+a,q+6 


=  G. 


m,n,pq 


-30- 


substitutlng  (3.10)  Into  (3.8)  and  (3.9)  and  rearranging,  the  final 
form  of  the  equations  Is  obtained  as: 
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The  calculation  of  the  coefficients  G  _  .is  discussed  in 

a  ,e  >Y><S 

Appendix  C.  The  final  step  is  to  solve  (3.10)  and  (3.11)  by  Gaussian 
elimination  and  backwards  substitution. 

One  point  deserve  further  comment.  Equation  (3.4)  was 
substituted  into  (3.3)  only  after  it  was  converted  into  the  form 
(3.10)  by  using  (2.22).  This  way  the  continuity  was  preserved, 
but  the  number  of  coupled  equations  was  reduced  ..proximately  by  one 
third.  As  a  check,  numerical  results  of  (3.10)  and  (3.11)  were 
compared  with  the  currents  calculated  from  the  corresponding 
equations  with  (3.10)  kept  as  a  separate  equation.  The  results  were 
identical  to  five  significant  decimal  figures. 

Having  determined  the  induced  currents,  it  is  now  possible 
to  calculate  the  far  field. 

3.4  The  Scattered  Field  and  Radar  Cross  Section 

The  far  field  is  derived  from  the  currents  using  the  Franz 
representation  and  the  far  field  approximation: 
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Frora  formulas  (3.5)  and  (3.6),  S  can  be  expressed  as: 


cos2  8$  sin  $s(Fy  cos  4>$  -  Fx  sin  <t>$) 
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e$  is  measured  from  grazing  and 


A  A  A  A 

r  *  x  cos  e$  cos  <j>s  +  y  cos  0S  sin  <f>s  +  z  sin  0S 

is  a  unit  vector  in  the  direction  of  the  observation  point 
(Fig.  1.1). 
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A  normalized  quantity  corresponding  to  the  scattered  far  field 
Is  the  bistatic  scattering  cross  section,  defined  as  (Crispin  and 
Siegel,  1968): 


a(r) 


11m 

r  -*  °o 


(3.17) 


Einc  Is  measured  on  the  surface  S.  For  a  plane  wave  Incidence 
(3.17)  can  be  written  as 


^  ~  |S(r) |2  .  (3.18) 

X2 

A  A 

The  value  of  a  corresponding  to  r  =  -k  is  called  the  backscattered, 
or  radar  cross  section. 

Equation  (3.17)  or  (3.18)  does  not  reveal  explicitly  the 
dependence  of  the  scattered  field  and  Its  polarization  on  the 
polarization  of  the  incident  field.  There  is  no  uniform  convention 
as  to  how  to  describe  this  dependence  in  the  general  bistatic  situation 
for  an  arbitrarily  shaped  body.  For  a  plate,  It  Is  convenient  to 
decompose  the  fields  into  components  parallel  and  perpendicular  to 
the  z  =  0  plane  (Fig.  1.1),  and  to  treat  separately  the  cases  of 
E-  and  H-polarizations,  l.e.,  Einc  or  H1nc  parallel  to  the  plane  and 
perpendicular  to  the  direction  of  incidence.  A  general  Incidence 
Is  a  superposition  of  these  two  cases.  Further  simplification 
results  from  choice  of  the  x-axIs  parallel  to  a  (b),  the  direction  of 
polarization  of  the  Incident  electric  (magnetic)  field.  With  this 
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choice  of  the  coordinate  system,  the  left-hand  sides  of  (3.11)  and 
(3.12)  are  ikY  exp(ikyn  cos  9^)  and  0,  respectively.  Restricting 
the  point  of  observation  to  the  plane  of  incidence,  i.e., 

♦  ■  ^  *  w/2,  the  direct  polarized  component  of  the  scattered 
field  becomes  actually  parallel  to  the  incident  field  (of  a  given 
polarization).  An  explicit  expression  for  (3.18)  is  then 
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for  direct  polarization.  The  cross  polarized  current  is  directed 
parallel  to 


u  * 


r  x  x  *  y  sin  9S  -  z  cos  es 


and 


CT 

A2 


Ic  12  p 

—  |Fy|» 


(3.20) 


(3.21) 


The  corresponding  integrals  are  given  by 

exp(-iky„  cos  ec)  (3.22) 
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3.5  Computed  Results  and  Comparison  with  Experiment 

Computations  were  carried  out  for  rectangular  plates  illuminated 
by  an  E-  or  H-polarizaed  plane  wave,  incident  in  one  of  the  principal 
planes  of  the  plate,  i.e.,  the  principal  axes  of  the  plate  coincide 
with  the  x-  and  y-axes.  The  direction  of  observation  was  confined  to 
the  plane  of  incidence. 

Computed  and  available  measured  results  for  1.16a  x  0.85a 
perfectly  conducting  and  resistive  plates  are  compared  in  Figs.  3.2 
through  3.5.  The  resistive  plate  is  a  composite  of  a  800  n/  square 
resistive  sheet  sandwiched  between  two  3.3  x  10“3  A  thick  dielectric 
plates.  The  relative  dielectric  constant  of  these  plates  is 
er  =  2.53. 

Figure  3.2  displays  the  normalized  radar  cross  section  as  a 
function  of  the  angle  of  incidence  8^  for  E-polarization  with  the 
field  parallel  to  the  long  side  of  the  plate.  Figure  3.3  presents 
the  corresponding  results  for  H-polarization.  Figures  3.4  and  3.5 
give  the  results  for  E-  and  H-polarizations  with  the  fields  parallel 
to  the  short  side. 

The  experimental  results  measured  over  an  angular  range  of 
180°  were  folded  into  a  90°  display,  revealing  some  asymmetries  about 
normal  incidence.  Examination  of  the  levels  of  the  peaks  at  that 
incidence  at  different  polarizations  show  inconsistencies. 

Nevertheless,  since  this  was  the  only  available  set  of  measured  data, 
it  was  used  for  making  a  comparison.  With  the  exception  of  the 
results  for  the  resistive  plate  in  Fig.  5.5  and  the  deep  nulls  in 
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all  the  other  cases,  computed  results  are  within  2  dB  of  the 
experimental  data  and  usually  much  closer.  There  are  two  reasons 
to  question  the  validity  of  the  experimental  data  of  the  particular 
case  in  Fig.  3.5.  At  normal  incidence,  the  maximum  level  is  several 
dB  higher  than  the  corresponding  peaks  at  other  polarizations,  whereas 
the  computed  values  are  mutually  consistent.  The  measured  lobe  near 
grazing  is  so  untypical  for  H-polarization  (theoretically  the  plate 
should  be  invisible  to  H-polarization  at  grazing  incidence),  that 
one  should  suspect  a  stray  reflection  from  the  experimental  setup. 

As  to  the  other  discrepancies,  all  computational  methods  consistently 
predict  deeper  nulls  than  are  actually  measured.  This  is  the  effect 
of  noise  filling  up  the  nulls.  Otherwise,  the  agreement  between 
computed  and  measured  results  is  quite  remarkable.  It  is  even  more 
so  considering  the  fact  that  owing  to  the  zero  edge  conditions,  the 
numerical  procedure  uses  an  effectively  smaller  plate.  A  zero  current 
is  enforced  on  half  cells  along  the  two  edges  of  the  plate 
perpendicular  to  the  current.  For  a  plate  divided  into  14  x  10  cells 
it  means  a  reduction  factor  of  13/14  or  9/10,  depending  on  the 
direction  of  incidence  and  polarization.  At  normal  incidence  the 
radar  cross  section  is  roughly  proportional  to  the  square  of  the 
area.  Thus,  the  corresponding  reduction  factors  for  the  cross 
section  are  -0.65  dB  and  -0.92  dB,  respectively.  This  brings  the 
results  into  even  closer  agreement  with  experimental  data. 

The  effective  resistivity  of  the  composite  plate  was 


computed  as  (Senior,  1981): 
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Req  ■  R  -  2ikd[cos  e1  -  YRUp  -  1)]J  (3.24) 

where  d  Is  the  thickness  of  the  dielectric  plates.  The  imaginary  part 

of  R  is  quite  small  and  the  value  of  R„„  at  normal  incidence  was 
eq  eq 

used  for  all  angles.  There  are  other  uncertainties  associated  with  the 
values  of  ep:  its  imaginary  part,  the  thickness  d  and  possible  air 
gaps  between  the  dielectric  plates  and  the  resistive  sheet.  It  is 
therefore  surprising  that  such  a  simple  model  produced  such  good 
resul ts. 

At  E-polarization  and  grazing  incidence  the  dominant  induced 
currents  on  a  perfectly  conducting  plate  are  edge  currents  along 
the  leading  and  trailing  edges  of  the  plate  (Knott  et  al . ,  1971).  A 
possible  qualitative  model  for  the  direct  polarized  scattered  field  is 
an  endfire,  two  element  array  of  wire  antennas  along  these  edges. 

This  model  is  supported  by  the  data  for  bistatic  scattering  in  the 
principal  plane  for  grazing  incidence,  presented  in  Fig.  3.6. 

Results  for  45°  and  normal  incidences  with  the  same  polarization  are 
shown  in  Fig.  3.7  and  Fig.  3.8,  respectively.  In  both  cases,  maximum 
scattering  occurs  in  the  specular  directions. 

Typical  current  distributions  for  grazing  and  normal 
incidences  are  shown  in  Figs.  3.9  through  3.11.  The  large  front 
and  rear  edge  currents  on  the  metallic  plate  are  clearly  visible  in 
those  figures.  The  currents  on  the  resistive  plate  do  not  possess 
such  singularities  and  their  variation  in  the  y-direction  resembles 
that  of  the  physical  optics  currents  (except,  of  course,  for  an 
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Fig.  3.10  Current  Distribution  Along  AA ’ ;  Normal  Incidence. 
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attenuation  factor).  The  variation  in  the  x-direction  is  of  the  same 
form  in  almost  all  cases;  therefore  only  one  typical  distribution 
is  presented.  Two  transverse  cross  sections  are  shown  for  the  current 
on  the  perfectly  conducting  plate  and  one  for  the  resistive  plate 
(on  which  there  is  very  little  variation  from  one  location  to  another 
along  the  y-axis). 

An  important  effect  of  the  cross  polarized  currents  is  the 
coupling  they  provide  between  the  front  and  rear  edges.  This 
coupling  mechanism  was  discussed  by  Knott  et  al .  (ibid).  Figure 
3.12  displays  the  behavior  of  the  edge  currents,  at  grazing  incidence, 
as  a  function  of  distance  along  the  edge.  Two  prominent  features 
of  this  graph  are  the  side  edge  current  buildup  and  the  nonsinu- 
soidal  nature  of  the  rear  edge  current.  These  characteristics  have 
been  detected  experimentally  by  the  above-mentioned  authors.  On 
the  resistive  plate  there  is  no  such  buildup  of  the  side  edge  current, 
which  is  small  compared  with  the  currents  on  the  leading  and  trailing 
edges . 

Results  for  scattering  by  plates  which  are  small  compared 
to  the  wavelength  computed ^by  the  procedures  described  in  this 
chapter  are  presented  in  Chapter  V,  after  the  topic  of  low  frequency 
scattering  has  been  discussed  in  Chapter  IV. 
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CHAPTER  IV.  LOW  FREQUENCY  SCATTERING  BY  PLATES 


4.1  The  Rayleigh  Series 

The  advantage  of  using  asymptotic  expansions  is  that  many 
analytic  and  semi-analytic  results  can  be  obtained  for  scatterers  of 
arbitrary  shape.  The  Rayleigh  series  is  such  an  expansion,  valid 
in  the  low  frequency  regime.  The  expansion  parameter  is  kD,  where 
D  is  a  characteristic  dimension  of  the  body.  The  series 
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and  the  corresponding  expansion  for  H  are  convergent  in  the  vicinity 
of  the  body  when  kD  <<  1 . 

The  individual  terms  of  the  Rayleigh  series  (4.1)  satisfy  the 
reduced  Maxwell's  equations 
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Stevenson  (1953)  has  offered  a  systematic  method  for  the 
determination  of  these  partial  near  fields.  The  essence  of  his 
method  is  that  the  terms  are  obtained  successively  as  solutions  of 
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potential  problems.  The  near  field  terras  are  then  used  to  produce  a 
multipole  expansion  for  the  far  field.  Later,  Kleinman  (1967)  showed 
that  in  some  cases  the  method  fails  to  produce  terms  beyond  the  first 
order  (ra  =  1 )  term  and  amended  the  procedure.  In  principle,  the 
revised  method  can  be  used  to  determine  the  partial  fields  of  any 
order.  However,  if  analytical  expressions  are  desired,  there  is  a 
difficulty  associated  with  the  resistive  (or  perfectly  conducting) 
boundary  condition.  For  every  ra,  is  the  sum  of  a  known  vector  and 
a  gradient  of  a  potential  function  (Kleinman,  1977).  The  boundary 
condition  yields  equations  for  the  tangential  derivatives  of  that 
potential.  In  some  cases,  such  as  m  =  0  or  bodies  of  special 
geometry,  these  equations  can  be  integrated  analytically  to  yield  a 
standard  Dirichlet  boundary  condition.  (Numerical  integration 
involving  an  open  line  integral  is  always  possible.)  In  some  instances 
at  least,  it  is  possible  to  derive  analytic  expressions  for 
the  first  order  far  field,  in  terras  of  zeroth  order  potentials, 
without  finding  the  corresponding  near  field  first. 

4.2  The  Scattered  Far  Field 

When  determining  the  scattered  far  field  at  low  frequencies 
the  Franz  representation  displays  a  clear-cut  advantage  over  that  of 
Stratton.  Franz's  formulation  is  the  "natural"  representation  for  the 
scattered  far  field,  since  that  field  is  generated  by  the  Induced 
currents  only.  The  induced  charges  have  no  direct  effect  on  the 
radiation  field.  One  near  field  term  is  needed  to  determine  the 
corresponding  far  field  term  if  the  Franz  representation  is  used, 
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but  two  near  field  terms  are  required  when  the  Stratton-Chu 
representation  is  employed  (Kleinman,  1967b). 

Using  the  Franz  formulation  with  electric  currents  only  and 
the  far  field  approximation,  the  scattered  field  is  given  by 
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Substituting  expansion  (4.1)  for  H  and  collecting  equal  powers  of  ik, 
the  field  is  expressed  as: 
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Further  simplification  is  obtained  by  the  vector  identity  (Jones, 
1964,  p.  531). 


,V  dS ’ 


J (n*v' „ 


V)r '  dS ’ 


(4.3) 


i 


( 


An  alternative  and  more  symmetric  form  is 


ES 


=  -k2 


ikr 


4ur 


7  rArA(pQ  +  ikp^ ) 


+  ZrA (m  +  i 


ikm^ 


+  0(k“)  ,  (4.8) 


which  is  obtained  by  further  applications  of  (4.3)  to  V  =  H  (r*r‘) 

l 

and  to  V  =  Ho(r*r')2.  However,  (4.4)  proves  to  be  more  convenient 
for  actual  calculations.  The  quadrupole  moments  in  (4.8)  are  given 
by 


P1  =  £l 


J  (n*E^)r'  dS'  -  ^  J  (n*Ej  (r*r '  )r'  dS * 


and 


•  7  J  r'A  (n^Hj  )dS'  -  f  Jr' \  (n.H0)(r.r' '  )dS' 


.  (4.9) 


Once  the  induced  currents  and  charges  are  determined,  the 
far  field  can  be  calculated  by  (4.4)  or  (4.8). 

It  is  now  necessary  to  deal  with  the  particular  geometry  and 
constitution  of  the  scattering  body.  The  following  sections  will 
focus  on  perfectly  conducting  and  resistive  plates.  The  perfectly 
conducting  plate  is  the  more  basic  structure  and  is  discussed  first. 


4.3  Perfectly  Conducting  Platelets 

The  first  two  terms  of  the  expansion  (4.1)  of  the  incident 
fields  (1.1)  and  (1.2)  are: 


|inc 
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Ho 
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The  currents  and  charges  on  a  flat  Infinitesimally  thin  plate  are 
given  as  differences  across  the  plate 

3  =  u\: 

and 

£■  =  z  -  E  | +  . 

£ 

Dipole  Moments.  The  static  fields  E^  and  R^  are 
Irrotational  (see  4.2)  and  can  be  obtained  from  potentials.  In 
the  case  of  E^.  E®  *  -v<j>q  where  <j>o  is  an  exterior  potential.  The 

A 

corresponding  incident  potential  is  -a*r,  and  from  the  boundary 
condition  =  -z.E^nc  on  the  plate  it  follows  that 

<j>0  =  a*r  +  c  on  S  ,  (4.14) 

where  c  is  a  constant  chosen  to  satisfy  the  zero  induced  charge 
condition 


(4.12) 

(4.13) 


/ 


!!oi+ 

3Z 


dS 


=  0 


(4.15) 


Because  of  the  form  of  (4.14)  it  is  convenient  to  write 


o 


(4.16) 


where  the  x4,  1*1 ,2,3,  are  Cartesian  coordinates  with  x  ■  2.  The 

1  3 

*  3 

potential  can  be  represented  as 

4  r  + 

♦o(f>  -  -Jw- 1  So  ds'  • 
s 


and  by  applying  the  boundary  condition  (4.14)  we  obtain 


a*1  + 

JT  Go  dS'  •  <4-17) 

s 

34>q 

which  is  an  integral  equation  from  which  to  determine  -r~r 

9  Z 

On  substituting  (4.16)  into  (4.5)  we  have 


p  *  cP*a  (4.18) 

0 

(Keller  et  al . ,  1972),  where  P  is  the  electric  polarizability  tensor 
with  elements 

Xj  dS  ,  (4.19) 

which  are  functions  only  of  the  geometry  of  S.  It  is  evident  that 
P^j  =  0  if  i,j  =  3,  and  from  (4.15),  (4.17'  and  (4.19) 

f  f  3$o|+  3*o  + 

P1j  =  JJGo  3#|_  aT_dSdS'  •  (4-20) 

The  tensor  is  therefore  symmetric,  i.e.,  P^j  *  P^. ,  and  has  at  most 
three  independent  elements. 
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The  procedure  for  the  magnetic  field  Is  similar:  H*  ■  -YVt|>o 
where  i/Q  Is  again  an  exterior  potential  In  terms  of  which 


7  *.i: 

~  r* 


(4.21) 


(Klelnman,  1973).  The  boundary  condition  on  the  plate  Is  z*Hc 
-z*H^nc ,  Implying 


=  b*z  on  S 


(4.22) 


Therefore  ij>  Is  normalized  as 
0 


1 1>  *  b»z  . 

To  0 


(4.23) 


The  representation 


r  3G 

♦2<?>  ■  J  JT  *21. 


then  leads  to  the  integral  equation 


-i  -  —  f  g0*»,+ 

3Z2  Js  0  0  - 


(4.24) 


from  which  I +  can  be  found.  In  terms  of  this 


mQ  =  Y(b-z)z  J 
s 


(4.25) 


i 
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implying  (Keller  et  al.,  1972)  a  magnetic  polarizabil Ity  tensor  M 
with  just  the  single  non-zero  element 

M  »  -  f  *3|+  <*S  .  (4.26) 

33  J  0  - 

s 

The  solution  of  equations  (4.17)  and  (4.24)  and  the  computation  of 
the  dipole  moments  are  discussed  in  Chapter  V. 

4.4  Resistive  Platelets 

The  resistive  plate  is  described  mathematically  by  (1.3)  and 

(1.4).  At  k  »  0  however,  the  condition  y  =  yQ  must  be  explicitly 

applied.  The  plate  is  therefore  invisible  as  regards  the  static 

magnetic  field,  i.e.,  HQ  *  0  and  (1.4)  then  implies  that  for  all  R, 

£  is  the  same  as  if  R  =  0.  Thus,  for  a  resistive  plate,  =  0 

o  o 

and  the  electric  dipole  moment  pQ  is  the  same  as  for  perfectly 
conducting  plates. 

To  determine  how  the  low  frequency  scattered  field  is 
dependent  on  R,  it  is  necessary  to  examine  the  first  order  terms 
in  (4.4).  Since  Hq  =  0  the  third  terra  of  (4.7)  is  zero.  To  evaluate 

the  second  terra,  z.H  has  to  be  determined  first.  From  Maxwell's 

l 

equations  it  is  found  that 
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3_ 

3Z 


Cz.H  ) 

l 


+  YzJ2^q) 


and  by  application  of  (1.4)  on  S 


Hence, 


on  S 


In  addition  =  0. 

The  Cartesian  components  of  are  themselves  exterior 
potentials  and  therefore 


H 

1Z 


G  dS  ’ 
o 


on  S 


A1  so. 


(4.27) 


(4.28) 


and  hence,  on  the  plate, 
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|r  (z.H  )  =  -—  /  z.H  |+  GQ  dS*  .  (4.29) 

32  1  3z  2  J  1  ’  0 

Using  (4.27),  (4.28)  and  the  fact  that  z.E^  =  -z.a  on  S,  the  left-hand 
side  of  (4.29)  can  be  evaluated,  and  the  solution  of  (4.29)  can  then 
be  written  as 


z.H  |+  =  -Yz. (z.aH3  [ +  +  1  b*z  u  (4.30) 

1  "  0  _  K 

where  u  is  a  tangential  vector  satisfying  the  integral  equation 


f 


6  d  V 
o 


(4.31) 


with  the  line  integral  taken  around  the  perimeter  C  of  the  plate, 
and  1 +  satisfies 

-1  -  J  *|l!  G0  «'  .  (4.32) 

3Z  s 

is  proportional  to  the  magnetostatic  potential  for  a  perfectly 
conducting  plate.  The  numerical  solution  of  (4.31)  requires  only 
a  simple  modification  of  the  program  used  to  solve  (4.32). 


1 
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Knowing  _  and  Q,  the  third  term  of  (4.7)  is 


f  z.H  |+  r*r'  dS ’  -  -z_(z..a)  f  <|>3|+  r*r'  dS ’ 

Js  1  '  JS  0  " 


+  Z 


H 


u(r')r*r'  dS '  .  (4.33) 


The  next  step  is  to  compute  the  first  term  on  the  right-hand 

side  of  (4.7).  Since  H  =  0,  E  is  the  gradient  of  an  exterior 

o  i 

potential  (Kleinraan,  1967),  i.e.,  E^  =  -v$  with 


r  !!i 

J  32 


dS  =  0  . 


Hence 


/_  /*  3<ji  + 

i-i}:  r-dS'  *  -  £  i,  J  ^  (xj  *  c,)  dS' 


1=1  S 


and  from  (4.14),  (4.16)  and  the  reciprocity  theorem  for  exterior 
potentials, 

2 

J  z-Ej*?'  dS'  =  -  £  xf  f  dS-  .  (4.34) 

s  1  1=1 
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The  zeroth  order  potential  t>Q  is  associated  with  an  Incident 
field  having  a  *  .  If  the  corresponding  first  order  scattered 

magnetic  field  is  fij.  Maxwell's  equations  imply 


z*  v*o  I  -  *  Zz*7-'^l-  » 


(4.35) 


and  when  this  is  substituted  into  (4.34),  the  vector  relation  in 
Jones  (1964,  p.  531)  can  be  used  to  give 

2 


I 


z-E.fcr'  dS'  =  Z  £  Xi  J  H1 1 +  dS*  . 


l  l 


i=i 


From  the  boundary  condition  (1.4) 


Lv<j>i  *  LlJnC  -  Rz^(z.Hi  J*) 


=  z 


~[a(k  r)  -  a  |  .  u(b-z)J 


and  hence 


i 

J  z* tj*  r  ‘  dS'  =  -Z  £x,  j  {(k*r' )a- yi||3  j*  a  -  ( b- z)u> 


*  -i  i  + 

z.H1  dS' 
l ' " 


(4.36) 
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Final  ly,  from  (4.30) 


z.fiV 


-Yz,.(z„xi}<j,3|*  +  £  z*(k.x.)u 


and  when  the  vectors  are  combined  to  eliminate  the  summation, 
f z*Ej+r'  dS1  =  f  /  ^ { k- r ' )z. (z.a)  -  (b*z)z. (z„u) 

s  Js  l 


+  (a-u)(z  k)  J  «|/3|^  -  |  z.(z_a)^[  +  ^3|  + 


-  ^  (z-k)  [  (k-r')a  u  -  (b*z)u*uj^ 


dS '  . 


(4.37) 


The  complete  first  order  contribution  is  therefore 
*  J  ^j(k  -  r)  r'  z„(z„a)  -  (b*z)z,(z.u) 

+  (a-u)(z„k)  j  z„  (z.a)  ^1* 

+  |  [ (r*r ' )(b‘Z)u  -  (z„k)  |  ( k* r ' )a* u  -  (b*z)u»uj 


(4.38) 


It  is  evident  that  F  has  components  only  parallel  to  the  plate. 
As  expected,  F  is  a  function  of  R,  and  the  presence  of  terms 
inversely  proportional  to  R  makes  explicit  the  discontinuous 
behavior  as  R  -*■  0.  For  any  given  R  0,  there  is  a  frequency  (which 


i 


( 
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depends  on  the  plate  dimensions)  below  which  F  no  longer  represents 
a  valid  correction  to  the  zeroth  order  contribution. 

A  case  of  some  interest  is  that  in  which  the  incident  magnetic 
vector  is  parallel  to  the  plate,  so  that  b*z  =  0.  Then,  z„k  = 

-(a-z)b  and 

F(r)  =  z.(z^a)  J  [  (k  -  r)-r*  -  j  dS' 

S 

+  ^  b(a-z)  J  £  k-r'_  R.  j  a-u  dS *  .  (4.39) 

JS  1 

If,  in  addition,  a*z  =  0  implying  normal  incidence  on  the  plate,  the 
second  integral  in  (4.39)  vanishes.  The  first  two  order  terms  in  the 
low  frequency  expansion  then  go  over  smoothly  into  the  result  for 
a  perfectly  conducting  plate  as  R  -►  0.  On  the  other  hand,  if 

A  A  A 

k*z  =  0  corresponding  to  grazing  (edge-on)  incidence,  a  is  in 
the  z  direction  (and  perpendicular  to  u)  and  F(r)  =  0.  This  is 
consistent  with  the  fact  that  any  electric  current  sheet  is  invisible 
to  a  plane  wave  at  grazing  incidence  for  perpendicular  polarization. 


CHAPTER  V.  NUMERICAL  SOLUTION  OF  THE  STATIC  INTEGRAL  EQUATIONS 


5.1  The  Electrostatic  Potential 

The  electrostatic  potential  is  derived  from  the  integral 
equation  (4.17).  Although  the  equation  is  not  explicit  because  of 
the  unknown  constant  c. ,  the  form  of  (4.17)  suggests  the  following 
splitting 


>1 


+  c , 


^i(2) 

0 


(5 


where  ^  ^ 1  ^  and  <f> 1  2  ^ 
o  o 


satisfy  the  integral  equations: 


and 


x 


i 


-1  = 


_  Gq  dS1 
Gq  dS' 


(5 

(5 


respectively.  With  this  definition  and  the  total  induced  charge 
condition  (4.15),  c^  can  be  determined  as 


(5 


Integral  equations  (5.2)  and  (5.3)  differ  only  in  the 


excitation  term  and  therefore  can  be  treated  alike.  They  can  be 
summarized  by  the  single  equation 
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G  dS  ’ 
o 


(5.5) 


The  weakly  singular  kernel  GQ  does  not  impose  continuity 
restrictions  on  the  solutions,  and  the  simple  pulse  expansions 

M  N 

z  z< 

p=l  q=l 


r  *  v 

ii 

)  a„  rect 

■x  -  xp  • 

l 

rect 

r  y  -  y*  ‘ 

*  -  c _ _ 

pq 

AX 

L  j 

(5.6) 


are  adequate, where  <t>  stands  for  <j>^  ^ ,  <j>^  or  </2^.  Equation  (5.5) 


is  converted  into  a  set  of  algebraic  equations 


M  N 


z  z 


a  G 

pq  0, 


m,n,p,q 


1  £  m  M 
1  <  n  <  N 


(5.7) 


when  sampled  by  the  impulse  functions  Wmn  =  5 (x-xm)5 (y-yn ) . 

G  is  similar  to  G  of  Section  3.3,  except  that  it 

°ra,n ,p,q  m,n,p,q 

involves  the  static  Green's  function  GQ  instead  of  G.  Explicit 

expressions  for  G  are  given  in  Appendix  C. 

m,n,p,q 

The  amount  of  numerical  labor  can  be  considerably  reduced  if 
the  plate  is  symmetric  with  respect  to  one  or  two  axes.  If  r  is  the 
reflection  of  r  with  respect  to  the  x..  axis,  then 
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3$ 


i(i) 


3Z* 


(?) 


3(f> 


id) 


3z 


(?) 


and 


3<j> 


(2) 


3Z 


(?) 


3$ 


(2) 


3Z 


(?) 


as  is  evident  from  (5.2)  and  (5.3).  Then,  (5.4)  implies  that 
ci  =  0  for  plates  with  at  least  one  axis  of  symmetry.  Furthermore, 
one  component  of  (5.7)  can  be  written  in  the  form 


a  ( G 

DQ 

y  m,n,-p,q 


m,n,-p,q 


1  <  m  <  M/2 
1  <  n  <  N 


when  that  axis  of  symmetry  is  the  x-axis  and  the  cells  are 
numbered  accordingly. 

Similar  equations  hold  for  plates  that  are  symmetric  with 

respect  to  the  y  axis  or  with  respect  to  both  axes. 

Convergence  of  the  solution  of  (5.7)  is  enhanced  if  the  cell 

size  is  allowed  to  decrease  toward  the  edge  to  better  accommodate 
-1  /2 

a  variation  of  p  (Meixner,  1972),  where  p  is  the  distance  from 
the  edge.  However,  the  possibility  of  treating  both  electrostatic 
and  magnetostatic  potentials  by  the  same  program  makes  that  variation 
of  cell  size  undesirable. 

Once  Eqs.  (5.7)  are  solved,  it  is  possible  to  compute  the 


elements  of  the  electric  polarizability  tensor  using  (4.19) 
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together  with  (5.1)  and  (5.4).  Numerical  results  are  presented  in 
Section  5.3. 


5.2  The  Magnetostatic  Potential 

The  integral  equation  (4.24)  from  which  the  magnetostatic 
potential  is  determined  is  much  more  singular  than  (4.17).  Conse¬ 
quently,  the  more  elaborate  scheme  of  Section  2.2.2  must  be  used. 

The  first  step  is  to  replace  the  second  normal  derivative  in 
(4.24)  by  the  surface  Laplacian  (32/3x2)+(32/3y2) ,  using  the 
definition  of  the  Green's  function:  v2GQ  =  -6.  Then,  the  tangential 
derivatives  are  represented  by  (2.22),  and  ^3 | +  itself  is 
represented  by  a  pulse  expansion 

M  N 

♦:i: s  Z  Z  bPq  rect 

p=i  q=l 


|  X  -  X  1 

y  -  yn 

_ _E 

rect 

_ a 

AX 

L  J 

(5.8) 


The  resulting  algebraic  equations  are 


1 


M  N 


4  £  Z  6pq 

p=l  q=i 


m+l/2,n,p,q 


+  G 

°m-l/2 ,n,p,q 


2Go 

m,n,p,q 


+  G, 


Ay  L  m,n+l/2,p,q  m,n-l/2,p,q 


-  2G„ 


m,n,p,q 


1  <  m  <  M 

(5.9) 

1  <  n  <  N 
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Symmetry  considerations  similar  to  those  presented  In  the 
previous  section  apply  to  (5.9).  Mzz  is  computed  using  (4.26). 

5.3  Selected  Numerical  Results 

Several  factors  determine  the  efficiency  of  the  numerical 
solution  of  (5.7)  and  (5.9).  The  translational  symmetry  of  Gmnpq 
allows  a  very  considerable  saving  of  matrix  filling  time  if  equal 
size  cells  are  used.  Because  the  matrix  elements  of  (5.9)  are 
combinations  of  terms  similar  to  the  matrix  elements  of  (5.7),  it  is 
possible  to  incorporate  the  solution  of  both  equations  into  a  single 
program,  leading  to  a  further  reduction  in  filling  time. 

Several  versions  of  the  computer  program  were  written.  Programs 
EPLT  and  APLT  solve  both  (5.7)  and  (5.9)  and  employ  uniform  cell  size. 
EPLT  is  for  plates  that  are  symmetric  with  respect  to  both  the  x  and 
the  y  axes  and  APLT  is  for  asymmetric  plates.  Program  DPLT  solves 
(5.7)  for  symmetric  plates.  DPLT  allows  several  optional  variations 
of  the  cell  size.  It  seems  that  the  quadratic  taper: 


3n* 

N(N  +  1  )(2N  +  1 


(5.10) 


is  the  best  compromise  of  decreasing  the  cell  size  near  the  edge 
without  increasing  the  middle  cells  too  much.  N  is  the  total  number 
of  subdivisions  per  half  plate  along  a  given  axis  and  n  =  1,2,...,N 
is  a  running  index  of  the  subdivisions,  counting  from  the  edge 
towards  the  middle. 


-70- 


Results  of  a  convergence  test  for  P  (*P  )  for  a  square  plate 

xx  yy 

are  displayed  in  Fig.  5.1.  Curve  No.  1  presents  results  computed  by 
EPLT,  and  curve  No.  2  the  corresponding  results  of  DPLT.  The 
convergence  rate  is  compared  with  that  reported  by  Okon  and 
Harrington  (1980,  p.  76).  They  used  triangular  cells  with  area 
coordinates  but  also  tried  uniform  cells  (curve  No.  3)  and  tapered 
cell  sizes  (curve  No.  4).  Since  their  program  was  intended  to  be 
as  general  as  possible,  they  have  avoided  taking  advantage  of  the 
symmetry.  To  make  the  comparison  in  terms  of  sampling  points  on  the 
plate,  their  reported  number  of  subdivisions  was  divided  by  2.  For 
the  same  N  their  matrices  are  four  times  as  large.  Other  adaptations 
of  their  results  were  multiplication  by  a  factor  2  and  an  interchange 
of  the  roles  of  the  electric  and  magnetic  dipole  moments,  since  they 
were  concerned  with  apertures. 

The  vertical  scale  of  Fig.  5.1  was  greatly  magnified  to 
emphasize  the  differences  in  the  convergence  rates  of  the  various 
methods.  As  expected,  the  methods  using  tapered  cell  size  converge 
much  faster  than  those  employing  uniform  cells.  The  quadratic  <.aper 
(5.10)  is  slightly  more  efficient  than  the  linear  taper  used  by 
Okon  and  Harrington  (ibid).  Their  results  computed  with  uniform 
cells  have  not  converged  yet,  but  seem  to  follow  the  same  curve  as 
the  results  computed  by  EPLT. 

There  was  a  question  as  to  how  well  our  programs  would  perform 
when  applied  to  geometries  that  do  not  conform  to  rectangular  cells. 
The  circular  disc  for  which  an  analytical  solution  is  available  was 
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used  to  test  the  performance.  Only  cells  whose  center  lay  inside  the 

circle  were  retained  and  the  effective  shape  was  polygonal  with 

serrated  edges.  Since  the  potentials  grow  indefinitely  towards  the 

edge,  this  procedure  could  conceivably  produce  large  errors. 

Normalized  values  of  P  (=P  )  obtained  by  EPLT  are  listed  in 

xx  yy 

Table  5.1  along  with  the  corresponding  errors.  Results  of  Okon 
and  Harrington  (ibid)  are  listed  as  well  for  comparison.  Although 
their  grid  matches  the  circular  shape  much  more  closely,  their 
results  are  better  only  when  tapered  cells  are  used.  Comparison 
of  the  results  of  EPLT  for  a  rectangle  and  a  disc  shows  that  the 
accuracy  was  reduced  from  a  small  fraction  of  a  percent  to  3  to  4  percent 
when  the  same  number  of  sampling  points  were  used. 

Figure  5.2  shows  the  convergence  of  Mzz  for  a  square  plate. 

The  curve  marked  by  6  =  0.5ax  was  produced  by  EPLT.  The  one  marked 
6  =  0.35ax  depicts  the  results  of  a  similar  program  where  the  nodes 
were  allowed  to  move  away  from  the  edge  of  the  cell,  as  shown  in 
the  insert.  A  third  line,  marked  6=0,  shows  the  results  of 
evaluating  the  second  normal  derivative  in  (4.24)  analytically  under 
the  integral  sign,  expanding  by  pulses  and  point  matching.  The 

rapidly  decreasing  rate  of  convergence  is  the  result  of  the  loss  of 
continuity.  The  utmost  importance  of  C°[S]  continuity  is  highlighted 
in  a  different  way  in  Fig.  5.3.  Mzz  is  plotted  as  a  function  of  6 
for  a  fixed  N  (=2).  Compared  to  the  value  displayed  in  Fig.  5.2, 
the  error  in  the  computed  value  of  Mzz  has  a  sharp  minimum  at 
6  =  0.5&X,  where  the  numerical  representation  is  continuous.  Most 


Okon  and  Harrington 
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Table  5.1:  Normalized  Dipole  Moments  for  the  Disc 


Method 


N 


mz2/a 


3/2 


Error 

(%) 


Pxx'A 


3/2 


Error 

(%) 


Error  (%) 


7 

9 

10 

12 


0.4604 

0.4614 

0.4608 

0.4620 


3.9 

3.7 

3.8 
3.5 


0.9082 

0.9184 

0.9249 

0.9297 


5.2 

4.1 

3.4 

3.0 


1.3 

0.6 

0.6 

1.0 


II 


Q) 

•O  I— 
•<-  «1 

(O 
> 


0.4506 

0.4572 


5.9 

4.5 


0.8904 

0.9038 


0.4754 

0.4766 


0.7 

0.46 


0.9644 

0.9570 


7.02 

5.64 


0.69 

0.09 


0.64 

0.41 


1  .35 
0.83 


Exact 


0.4788993 


0.9577986 
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notlceable  Is  the  jump  occurring  between  &  ~  0.4995ax  and  6  *  0.50ax. 
The  same  behavior  was  observed  for  larger  values  of  N,  where  the  scale 
was  compressed  due  to  the  convergence.  The  convergence  rate  with 
6  3  0.5ax  Is  comparable  to  that  of  Okon  and  Harrington  (ibid)  with 
variable  cell  size.  The  opposite  direction  of  convergence  is 
attributed  to  their  use  of  different  basis  functions  which  under¬ 
estimate  the  potentials  over  most  of  each  cell. 

The  accuracy  of  the  results  for  other  geometries  was  again 
tested  using  a  disc,  and  the  results  are  summarized  in  Table  5.1. 

The  accuracy  is  comparable  to  that  of  Okon  and  Harrington  (ibid)  with 
uniform  cells.  The  accuracy,  however,  depends  on  the  particular 
choice  of  the  grid  and  care  must  be  exercised  in  selecting  this 
grid.  Meshes  that  consistently  overestimate  or  underestimate  the 
area  of  the  plate  should  be  avoided. 

Numerical  results  for  the  dipole  moments  of  several  plates 
are  presented  in  Table  5.2.  The  shapes  of  these  plates  ar~  displayed 
in  Fig.  5.4.  For  each  shape  the  normalized  dipole  moments  are  given 
for  several  values  of  the  length  to  width  ratio  (L/W).  For  a  given 
value  of  L/W,  the  dipole  moments  of  plates  that  are  convex  in  the 
plane  show  very  little  sensitivity  to  the  shape.  Shapes  that  are 
concave  in  the  plane  deviate  from  this  rule,  but  the  largest  deviation 
is  associated  with  doubly  connected  shapes.  These  observations 
are  supported  by  the  data  of  Oe  Smedt  and  Van  Blade!  (1980).  For  such 
shapes,  normalization  to  (area)3/2  is  no  longer  appropriate  and. 


in  fact,  the  polarizabil ity  tensor  elements  diverge  as  the  ratio 


rectangl e 


rhombus 


Fig.  5.4  The  Various  Shapes  Used  in  Table  5.2 
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Shape 

L/W 

rectangle 

1 

2 

5 

rhombus 

2 

5 

cross 

1 

2 

5 

frame 

1 

2 

5 

Table  5.2:  Dipole  Moments  of  Various  Plates 


n 


Mzz/A 


0.4479 

0.4474 


0.4165 

0.4148 

0.4144 


0.3129 

0.3117 

0.3113 


0.4050 

0.4052 


0.3028 

0.3027 


0.4290 

0.4268 


0.3007 

0.2994 


0.2490 

0.2445 


0.2423 


0.2201 


0.4027 

0.3981 


P  V/AJ 

xx 


0.9889 

0.9972 


.6 
.6 
0.6185 


0.3513 

0.3566 

0.3584 


0.5788 

0.5841 


0.3289 

0.3317 


0.9816 

0.9966 


0.5954 

0.6045 


0.3393 

0.3442 


0.9477 


0.5490 


1 .095 
1.116 


0.3040 

0.3070 
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of  the  area  of  the  plate  to  the  area  of  a  singly  connected  plate  of 
the  same  exterior  contour  becomes  smaller. 

In  the  case  of  a  solid  body  a  variety  of  isoperimetric  and 
other  bounds  on  the  tensor  elements  have  been  developed  (Kleinman 
and  Senior,  1972),  and  may  serve  to  adequately  approximate  the 
elements.  For  a  plate  the  more  useful  estimates  are  empirical  in 
nature  (De  Smedt  and  Van  Blade!,  1980),  but  bounds  can  be 
rigorously  established  from  the  fact  that  for  a  solid  body 

3  3  3 

P  1  Z  I  Pij(a-ii)(a-ij)  and  M  =  £ 

i=i  j=l  i=i  j=1 

are  increasing  set  functions  (Schiffer  and  Szego,  1949).  By 
regarding  the  plate  as  the  limit  of  a  right  prism,  we  have 

P'  .  <  P  .  <  P" 

u  —  1 1  -  i  i 

(1-1,2)  , 

M'  <  M  <  M" 

3  3  -  3  3  -  ?  3 

where  the  single  and  double  primed  quantities  refer  to  the  largest 
inscribed  and  smallest  circumscri bed  circular  disks  respectively. 
For  a  disk  of  radius  a,  P  =  P  =  16a3/3  and  M  =  8a 3 /3 .  These 

1  1  2  2  3  3 

bounds  are  not  applicable  for  doubly  connected  shapes. 

It  is  evident  from  Table  5.2  and  other  data  published  in  the 

literature  (De  Smedt  and  Van  Bladel,  1980;  Okon  and  Harrington, 

1980),  that  P  and  M  decrease  with  increasing  L/W,  whereas  P 

xx  ul.  yy 


V  M .  . (b*  x . ) ( b* x  . ) 

L-  i  j  i  j 


-80- 


increases  (y  is  the  long  axis  of  the  plate).  The  bounds  in  terras  of 
inscribed  and  circumscribed  circular  disks  become  useless  for  large 
values  of  L/W.  Other  bounds,  in  terms  of  circular  disks  of  the  same 
area,  can  be  conjectured,  but  no  mathematical  proof  is  yet  available 
for  them,  and  there  is  always  a  risk  that  some  shape  exists  that 
violates  them. 

As  an  additional  test,  a  comparison  was  made  of  the  dynamic 
and  the  static  programs.  It  is  possible  to  compute  (less  efficiently) 
the  dipole  moments  by  the  program  for  dynamic  scattering.  Making 
the  plate  very  small  compared  to  tne  wavelength  and  using  E-polarization 
at  grazing  incidence  so  that  b-z  =  a-x  =  1,  the  dipole  moment  m0  is 
given  by 


*0  *  J  f  r' J  dS-  =  \  f  (x'Jy  -  y'JjdS1  , 
s  J  s 


and  by  using  (3.5)  and  (3.6), 


M  N-i 
C  <; — > 


M-i  N 


'  ?aXay{£  £  X"Vn+,/2  ■  l 


zz 


5.10) 


Similarly, 


M  N 


xx 


*  /  f x' ds'  ■  ^  Z  Z  t1  x» 


m=i  n=l 


iZ 


ixAy 


M-:  N 

I  ZJ- 

m  =  1  n  =  ’. 


'5.11  ) 


m ,  n + 1  /  2 


F/B  20/14 
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Computed  results  for  square  plates  are  shown  In  Fig.  5.5. 

The  values  of  P  (*P  )  for  both  perfectly  conducting  and  resistive 

aa  yy 

plates  approach  the  static  limit  as  U  decreases.  The  magnetic  dipole 
moment  of  the  metallic  plate  Is  close  to  the  static  limit  and  the 
magnetic  moment  of  the  resistive  plate  disappears.  These  results 
were  predicted  In  Chapter  IV,  and  are  used  as  an  additional 
verification  of  the  computer  program. 
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CHAPTER  VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  problem  of  electromagnetic  scattering  by  resistive  and 
perfectly  conducting  rectangular  plates  was  solved  numerically.  The 
numerical  techniques  developed  for  this  problem  can  be  easily 
applied  to  other  problems  involving  E-fleld  Integral  equations, 
e.g.,  wire  antennas,  strips  and  non-pi anar  open  shells.  Impedance 
plates  can  be  treated  as  well  by  superimposing  resistive  and  conductive 
plates. 

The  numerical  procedure  was  tested  by  comparison  with  measured 
data  and  with  analytical  and  computed  results  for  very  small  plates. 

The  accuracy  of  computed  radar  cross  sections  for  plates  of  about 
a  square  wavelength  In  area  Is  1  to  2  dB.  The  qualitative  behavior 
of  the  Induced  currents,  especially  the  edge  currents,  checks  well 
with  currents  measured  on  slightly  smaller  plates. 

Numerical  representations  were  developed  for  a  variety  of 
operators.  Comparison  of  Eqs.  (3.11)  and  (3.12)  with  the  Cartesian 
components  of  Eq.  (3.1)  reveals,  a  posteriori,  the  following  numerical 
representations  for  the  differential  operators  a2 /3x3y  and 
32/ayax: 


lit. 

axay 


-  f<V*n>  *  ,(V»n  ' 


(6.1) 
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and 


>iq  .  _ 

»H  (*..*„>  A“y 


fe  ['(V>„^)  -'(*■-  *•' n*¥) 

-  f<V*n>  +  f(xm  ’2**  *n> 


(6. 


that  preserve  the  continuity  of  f  given  by  (2.22)  and  (2.23),  and 
with  these  representations  It  Is  possible  to  deal  with  E-fleld 
Integral  equations  derived  from  either  the  Franz  or  Stratton-Chu 
formulations. 

At  present,  the  computer  program  has  two  major  limitations: 
It  handles  only  rectangular  plates  and  the  maximum  area  of  these 
plates  Is  restricted  to  about  a  square  wavelength.  Both  short¬ 
comings  can  be  removed  or  alleviated.  There  are  several  possible 
courses  of  action  which  enable  other  shapes  to  be  handled.  The 
most  straightforward  one  Is  to  try  to  approximate  that  shape 
using  rectangular  cells.  This  method  was  used  In  the  static 
program  and  produced  results  for  the  dipole  moments  to  an  accuracy 
of  a  few  percent.  The  radar  cross  section  Is  a  weighted  Integral 
of  the  Induced  currents  and  Is  thus  expected  to  be  somewhat 
stationary. 

Another  approach  Is  to  map  the  rectangular  mesh  onto  a  grid 
with  general  quadrilateral  cells,  to  fit  more  accurately  the  given 
shape.  Experience  with  the  static  programs  may  discourage  one  from 
taking  this  course.  The  cost  of  running  the  program  that  employs 


a a 


-85- 


unequal  cells,  was  higher  than  that  of  the  other  program,  using 
equal  cells,  with  matrices  about  twice  as  large, (144  x  144) 
compared  with  (64  x  64).  The  reason  is  that  the  latter  program 
takes  advantage  of  the  translational  symmetry  of  the  matrix  elements 
to  reduce  the  number  of  different  elements  to  M  x  N,  where  N  and  N 
are  the  number  of  subdivisions  along  the  x  and  y  axes,  respectively. 
This  cannot  be  done  with  unequal  cells  and  the  number  of  matrix 
elements  to  be  computed  rises  to  (M  x  N)2.  The  fill  time  of  the 
corresponding  matrix  exceeds  the  Inversion  time  of  the  much  larger 
matrix  In  the  other  program. 

A  possible  compromise  Is  to  use  special  shapes  for  edge 
cells  alone.  Since  the  number  of  edge  cells  Is  but  a  fraction 
of  the  total  number,  the  matrix  fill  time  will  increase  less  sharply. 
Another  feature  that  can  be  Incorporated  Into  a  special  treatment 
of  the  edge  cells  Is  to  use  singular  elements  (Zlenklewlcz,  1977, 
p.  670;  Jones,  1979,  p.  506).  This  may  somewhat  Increase  the  rate 
of  convergence,  and  Its  effect  will  be  more  pronounced  In  the  static 
programs,  where  the  number  of  required  sampling  points  does  not 
depend  on  the  size  of  the  plate. 

Non>planar  open  shells  can  be  treated  by  mapping  the 
rectangular  grid  onto  a  curved  mesh.  A  simple  example  Is  Gllsson 
and  Wilton's  (1980)  treatment  of  a  bent  plate.  In  general,  such  a 
mapping  will  result  In  unequal  cells. 

Unlike  the  shape  limitation  that  can  be  relatively  easily 
removed  by  a  modification  of  the  existing  program,  an  Increase 
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In  the  maximum  plate  size  that  can  be  handled  requires  Improvement 
In  the  matrix  Inversion  techniques.  The  Interaction  matrix  for 
the  coupled  Integral  equations  Is  a  dense  matrix.  When  equal  cells 
are  used,  the  matrix  Is  symmetric  and  there  Is  an  additional  hidden 
symmetry  because  the  kernel  of  the  equation  Is  of  a  convolution 
type.  The  Interaction  matrices  of  one-dimensional  structures  with 
convolution  kernels  are  of  Toeplltz  form.  However,  there  Is  no 
way  In  which  the  corresponding  Interaction  matrix  of  a  two-dimensional 
structure  can  be  put  Into  that  form.  Some  progress  In  Inversion 
techniques  for  matrices  that  are  close  to  Toeplltz  matrices  has 
been  reported  by  Frledlander  et  al .  (1979),  but  the  technique  Is 
not  applicable  to  plates  with  variable  resistivity  because  the 
elements  on  the  main  diagonal  of  the  corresponding  Interaction 
matrix  are  all  different.  Matrices  of  this  type  may  have  to  wait 
for  hardware  Improvements  that  will  permit  the  handling  of  much 
larger  matrices. 
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APPENDIX  A.  DERIVATION  OF  THE  REPRESENTATIONS  FROM  THE 
DEFINITIONS  OF  THE  DYADIC  GREEN'S  FUNCTION 


The  electric  and  magnetic  dyadic  Green's  functions  are 
defined  (Tal,  1973)  by: 


’•8e  ’  \  <*•’> 

and 

=  k2§e  +  !<s(r  -  r')  .  (A. 2) 


They  also  satisfy 

-  i-v[l«(r  -  ?’)]  (A. 3) 

k2 

7-<L  *  0  .  (A. 4) 

ra 

Equations  (A.l)  and  (A. 2)  can  be  replaced  by  the  uncoupled. 


and 


second  order,  vector  Helmholtz  equations 
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Levlne  and  Schwinger  (1950)  concluded  from  (A. 3),  (A. 4), 
(A. 5),  (A. 6)  and  the  Helmholtz  equation  for  the  scalar  Green's 
function 

V2G  +  k2G  *  -S(r  -  r') 


that 


I  + 


~ —  77 1 

k2  J 


G 


(A.  7) 


and 

Gm  -  ^(fG)  .  (A. 8) 


Equations  (A. 5)  and  (A. 6)  are  functional  equations  In  a 
distribution  space  (and  so  are  (A.1)-(A.4)).  Before  interpreting 
them  In  that  way.  It  Is  necessary  to  specify  the  distribution  space. 
Gelfand  and  Shilov  (1964)  noted,  in  the  preface  to  the  second 
volume  of  their  book,  that  different  distribution  spaces  should  be 
used  for  solving  different  types  of  problems.  It  Is  equally  true 
when  solving  scattering  problems  with  different  geometries.  For 
scattering  by  finite  smooth  bodies,  the  most  familiar  distribution 
space  can  be  used.  This  Is  the  space  of  linear  functionals  over 
the  space  of  infinitely  differentiable  vector  functions  3(r'),  with 
finite  support  on  the  scatterer.  A  different  space  Is  required  for 
finite  scatterers  with  edges.  The  test  functions  are  differentiable 
and  bounded  only  In  the  Interior  of  the  support,  and  some  of  their 
components  may  diverge  towards  the  edges  as  p_t  (0  <  t  <  1) 

(Melxner,  1972).  A  countable  set  of  norms  for  this  space  is 
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1(3|(  ■  sup  J  pp3(r* )  |  p  *  1,2 . 

r'  £  S 

With  these  norms  the  proof  that  the  test  function  space  Is  a 
fundamental  one  and  the  construction  of  the  dual  space,  are 
exactly  analogous  to  those  of  Gelfand  and  Shilov  In  Vol .  2,  Ch.  II 
and  Vol.  3,  Ch.  I. 

Other  distribution  spaces  should  be  used  for  Infinite  and 
semi -Infinite  bodies  with  or  without  edges,  but  they  will  not  be 
discussed  here. 

Equation  (A. 5)  should  be  understood  as  a  functional  over 
the  test  functions  J(r'): 

<v*vAGe,3>  -  k2  <Ge,J>  *  <16 ,3>  . 

Since  G  Is  a  convolution  type  distribution,  the  differential 
e 

operator  can  be  brought  outside  to  give 

(vAVA  -  k2)  <GeJ>  »  3  *  -  j—  (va7a  -  k2 ) E  . 

The  last  equality  Is  obtained  from  the  vector  Helmholtz  equation 
for  E.  §e  satlslfes  the  radiation  condition  and  the  same  boundary 
conditions  as  E.  The  uniqueness  of  the  solution  of  the  Helmholtz 
equation  therefore  Implies 


J 
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E  -  IkZ  <Gfi,J>  -  ikZ 


[' +  7  +’ 


J>  . 


The  distribution  G  corresponds  to  the  ordinary  function  G(r.r'); 


hence 


i,3>  =  J  G(r,r' 


)J(F')  dS' 


As  far  as  the  scattered  field  alone  is  concerned,  the  body 
is  replaced  by  the  induced  currents,  and  the  free  space  Green's 
function  should  be  used.  Hence 


Es  =*  ikZ 


G  J  dS' 


(A. 9) 


and  a  similar  interpretation  of  (A. 6)  leads  to: 


Hs  =  VA  /  GJ  dS'  . 


(A. 10) 


Equations  (A. 9)  and  (A. 10)  are  identical  with  the  Franz 
representation  for  the  scattered  field  when  only  electric  currents 
are  present.  Magnetic  sources  can  be  treated  analogously  and  the 
complete  representation  is  then  obtained  by  superposition. 
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APPENDIX  B.  PROOF  OF  C°[fl]  CONTINUITY  OF  THE  SUBSECTIONAL  BASIS 

Using  the  notation  of  Fig.  B.l,  the  coefficients  a. 
expansion  (2.21)  for  cell  I  are  given  by: 


*  f CO ) 


cx=2 

s 


•  *.■  2 


a  *  — < —  [f(6)  -  f(8)  ♦  f(5)  -  f ( 7 )  +  2f(4)  -  2f(2)] 
7  Ax2Ay 


a  =  — ~  Cf(6)  -  f ( 8)  -  f{5)  +  f(7)  +2f(l)  -  2f(3)] 
8  AxAy2 


a  *  - -  <  f(5)  +  f ( 6 )  +  f ( 7 )  +  f(8)  -  2[f(l )  +  f(2) 

9  AxzAy2  l. 


+  f(3)  +  f  (4 )]  +  4f(0) 


(B.l) 
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* 

A  similar  expression  with  coefficients  and  variables 
*  * 

£  3£=x-xra>  n  *  y  -  yn+1  holds  for  cell  II.  Taking  the 

difference  a=  f_  _(e,n)  -  f_  n  )  along  the  common  line  of 

in  9  n  m,  n  ■  1 

the  two  cells  yields 

A  =  (a  -  a*)  ^  +  (a  -  a*)  +  (a  -  a*)  +  (a  +  a  *) 

1  1  6  66^  22  4  4  ^ 

2 

+  (a  -  a*)  +  (a  -  a*)  +  (a  +  a*)  4^  +  (a  -  a*) 

8  8  4  5  5  7  7  2  9  9'  4 

=  A  +  Be  +  C£2  .  (B.2) 


Using  (B.l )  and  the  corresponding  expressions  for  cell  No.  II  In 
(B.2),  the  coefficients  A,B  and  C  are  obtained: 
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Ax2-C  •  2£f(3)  +  f(l)  -  2f(0)]  -2[f(12)  +  f(10)  -  2f(9)J 

+  £f(6)  -  f(8)  +  f(5)  -  f(7)  +  2f(4)  -  2f(2)J  +  [f(14)  -  f(5) ' 

+  f(13)  -  f(6)  +  2f (2 )  -  2f(n )]  +  [f(5)  +  f(6)  +  f(7)  +  f(8)] 

-  2[f(l)  +  f(2)  +  f(3)  +  f(4)]  +  4f(0)  -  [f(13)  +  f(14)  +  f(6) 

+  f(5)]  +  2[f(10)  +  f(ll)  +  f( 1 2 )  +  f(2)]-  4f(9)  *  0  . 

N 

Therefore  f  ■  ^  fmn  Is  continuous  across  the  cell  boundary. 
N-i 


APPENDIX  C.  EVALUATION  OF  THE  MATRIX  ELEMENTS 


No  closed  form  Is  available  for  the  Integrals 


G 

u»v,tc, 


X 


X  +&X/ 2  y.+Ay/2 


f  f  6(xp.yv,x' .y'Jdx*  dy' 


Xk-AX/2  y*-Ay/2 


of  Section  3.3. 

When  the  observation  point  Is  not  In  the  Integration  cell, 
the  Integrand  Is  smooth  and  regular.  Numerical  Integration  by  2x2 
Gaussian  quadrature  is  optimal  (Zienklewlcz,  1977;  p.  280).  The 
self-cell  term  Is  approximated,  when  the  cell  Is  small  enough,  as 
follows : 


mnmn 


If 


Jk|r  -r1 
e  ran 


«-  .  i-  f  -l^'1 


S  lrmn-p'l 

mn 


ran 
+  IkAxAy 


dS’ 


(C.l) 


mnmn 


Gq  Is  the  self-cell  terra  of  the  mtarlx  of  the  electrostatic 

m,n,ra,n 

equation  (5.7).  This  Integral  can  be  evaluated,  using  cylindrical 
coordinates  and  the  notation  of  Fig.  Cl (a): 


J 


Fig.  C.1  Geometries  for  the  Calculation  of  the  Matrix  Elements 


»  4 


TB,n,m,n 


*o  Ax/2sec*  tr/2  Ay/2csc* 

r  d ♦  r  dp  +  f  <*♦  r 

K  J0  J*o  Jo 


•  2  ^  Ax  •  slnh-1  fj  ♦  Ay  •  csch-1  fj  |  • 


(C.2) 


The  other  cases  of  Fig.  C.l  are  evalauted  In  a  similar  fashion: 


2Ax  •  slnh- lM  ♦  Ay  csch-1 


am,n,m+1/2,n 


(C.3) 


as  Illustrated  In  Fig.  C.l(b),  and 


rn.n.K.n 


o  •  slnh-1  J  -  8”  slnh-1  J 


► 

i-1  g-  ♦  Y  csch-1  £  -  csch 
- 


C.4) 


where 


l*m 


x  + 

ic 


AX 

T  ’ 


l*m -  \>  -  r 


and 


Aj£ 

Y  *  2* 


as  Illustrated  In  Fig.  C.l(c)  and 
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ro,n,ic,X 


2  ^ajsin"1  -  slnh"1  ♦  B^slnh"1  -  slnh"1  ij 
♦  ^[csch*1  *•  -  csch"1  |-]  +  fi^csch"1  j  -  csch"1  (C.5) 


where  a  and  s  are  as  before  and 


and 


y  »  |yn  -  yx I  +& 

6  •  |yn  -  ^ 
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Coeiputed  results  for  square  plates  are  shown  in  Fig.  5.5. 

The  values  of  P  (»P  )  for  both  perfectly  conducting  and  resistive 

jm  jj 

plates  approach  the  static  Unit  as  W  decreases.  The  magnetic  dipole 
noaent  of  the  metallic  plate  Is  close  to  the  static  limit  and  the 
nagnetlc  moment  of  the  resistive  plate  disappears.  These  results 
were  predicted  In  Chapter  IV,  and  are  used  as  an  additional 
verification  of  the  computer  program. 


i 


t 
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